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Calculation of Thermally-Induced Displacements 
in Spherically Domed Ion Engine Grids 

George C. Soulas 

National Aeronautics and Space Administration 
Glenn Research Center 
Cleveland, Ohio 44135 

Abstract 

An analytical method for predicting the thermally-induced normal and tangential displacements of spherically 
domed ion optics grids under an axisymmetric thermal loading is presented. A fixed edge support that could be 
thermally expanded is used for this analysis. Equations for the displacements both normal and tangential to the 
surface of the spherical shell are derived. A simplified equation for the displacement at the center of the spherical 
dome is also derived. The effects of plate perforation on displacements and stresses are determined by modeling the 
perforated plate as an equivalent solid plate with modified, or effective, material properties. Analytical model results 
are compared to the results from a finite element model. For the solid shell, comparisons showed that the analytical 
model produces results that closely match the finite element model results. The simplified equation for the normal 
displacement of the spherical dome center is also found to accurately predict this displacement. For the perforated 
shells, the analytical solution and simplified equation produce accurate results for materials with low thermal 
expansion coefficients. 
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Non-homogeneous solution variable 
Non-homogeneous solution variable 
Constants for quadratic temperature equation 
Non-homogeneous solution variable 
Boundary condition constant 
Displacement boundary condition constant 
Elastic modulus 

Flomogeneous solution for displacement general solution 
Shell thickness 
Spherical dome height 
Defined variable 

Circumferential bending moment per unit length 
Meridional bending moment per unit length 
Series summation variable 
Defined variable 

Circumferential longitudinal force per unit length 
Meridional longitudinal force per unit length 
Legendre function 

Legendre, or spherical, function of the first kind 

Shear function general solution 

Shear function non-homogeneous solution 

Shear function homogeneous solution 

Shear force per unit length 

Spherical shell radius of curvature 

Chord radial location 

Spherical shell chord radius 

Shell change in temperature distribution 


NASA/TM— 2006-2 14046 


1 



V 
w 
W 0 
z 
a 
A 

Ar 0 

See 

S<pip 

F 0 

£ ee 

e ° 

c w 

0 

Ke 

K<p 

V 

Oee 

<V 

^(ptp 

i 

To 

9 

<Pe 

1 

X 

X 

V 


Shell tangential displacement 
Shell normal displacement 
Displacement w at (p = 0 

Radial direction in a spherical coordinate system 
Thermal expansion coefficient 
N on-homogeneous solution variable 
Spherical shell change in chord radius 
Circumferential strain 
Meridional strain 

Circumferential middle surface strain 
Meridional middle surface strain 
Circumferential direction 
Circumferential strain equation variable 
Meridional strain equation variable 
Poisson’s ratio 

Circumferential (or hoop) stress 
Shear stress 
Meridional stress 

Shell temperature distribution in meridional direction 
Shell temperature distribution through thickness 
Meridional direction 
Meridional spherical shell edge 
Displacement function general solution 
Displacement function non-homogeneous solution 
Displacement function homogeneous solution 
Boundary condition constant 


I. Introduction 

A key ion engine component is the acceleration system, or ion optics assembly. The ion optics assembly is 
typically made up of two perforated electrodes, labeled the screen and accelerator grids. Critical dimensions that 
directly affect ion optics’ performance and service life are the gap between the grids and the grid-to-grid aperture 
alignment. 

During ion thruster operation, the grids are heated by power deposition from the plasma and conducted and 
radiated heat from other thruster surfaces. This heating typically results in a temperature distribution along the grid 
surface that is higher at the center of the beam extraction area and cooler at the beam extraction edge. To preclude 
buckling and control the thermal expansion of large area grids, the thruster grids are typically manufactured with a 
spherically domed shape. This shape significantly reduces compressive stresses that can cause buckling and causes 
both grids to thermally expand in the same direction to provide some control over the grid gap. 

However, each grid is also differentially heated during ion thruster operation. The resulting uneven 
displacement of each grid during thruster operation causes the gap between the grids and the radial grid-to-grid 
aperture alignment to change from their preset values. Knowledge of the grid gap and aperture alignment during 
thruster operation is essential in properly assessing ion optics’ performance and service life. 

Direct measurements of grid gaps during thruster operation have been made, however, these measurements are 
difficult to do and consume considerable resources (refs. 1 to 5). The only known measurements of grid-to-grid 
aperture alignment during thruster operation are those of reference 3. Radial grid-to-grid aperture misalignment 
during thruster operation has typically been inferred by non-axisymmetric erosion of the accelerator apertures 
resulting from extended thruster operation. Accurately predicting the displacements of each grid both normal and 
tangential to the dome surface during thruster operation would be of value for predicting the grid gap and grid-to- 
grid aperture alignment, respectively. This is because grid gap and grid-to-grid aperture alignment changes at a 
particular location are merely the difference of the normal and tangential displacements, respectively, of each grid. 
Such a predictive technique would not only be of value for existing ion thruster designs, but also for future designs. 
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Predicting these displacements requires knowledge of two critical elements. The first element is the temperature 
distribution across each grid. The second element is the resulting, thermally-induced displacements of the grids from 
these temperature distributions. The latter element is the subject of this paper. 

Predicting grid displacement during thruster operation from a given temperature distribution has been the topic 
of several past investigations, only a few of which are included in the reference section of this paper (refs. 6 to 8). 
All of these studies, however, relied on finite element structural analysis computer codes, or finite element models, 
to determine the thermally-induced displacements of the grids. While such a technique can be effective in accurately 
determining these displacements, it lends no direct insight as to how grid design parameters affect grid 
displacements. Although these models can be used indirectly in modeling experiments where design parameters can 
be varied to determine their impact on displacements, such an approach is inefficient due to the broad parameter 
space. A more useful predictive tool would be an analytical method. 

This paper describes an analytical method for calculating the thermally-induced displacements of spherically 
domed grids. Displacements determined include those both normal and tangential to the grid surface for grid gap 
and grid-to-grid aperture alignment determinations, respectively. The theory of thin-walled spherical shells under 
thermal loads and the resulting stresses and moments are reviewed. The appropriate boundary conditions are applied 
for the determination of constants. Afterwards, normal and tangential displacements are derived. A simplified 
equation for the normal displacement at the center of the spherical dome is also derived. The process for including 
the effects of shell perforation is presented. Finally, the resulting analytical model is compared to the results from a 
finite element model. 


II. Background — Thermal Stresses in a Thin-walled Spherical Shell 

The theory of thin-walled shells of revolution under thermal loads varying both through the thickness and along 
surface was originally reported by Nowacki (ref. 9), and later by Johns (ref. 10). A summary of the theory and the 
analytical solutions for thermal stresses in thin-walled shells from reference 9 is presented below. It will serve as 
a background necessary for determining the thermally-induced displacements of a spherical dome and provide 
equations for determining thermally-induced stresses and moments. The notation used in reference 9 will be 
adopted here. 

In the following analysis, only spherical shells of uniform thickness under an axisymmetric thermal loading 
(i.e., constant along the circumference) will be considered. It is assumed that the wall thickness is very small in 
comparison with the radius of curvature and that the normal stress through the wall thickness is zero, both of which 
are commonly assumed in thin shell theory (ref. 11). The change in curvature resulting from the thermally-induced 
displacements is assumed to be negligibly small and that buckling does not occur. Furthermore, it is assumed that all 
deformation is fully elastic and obeys Hooke’s law. Finally, it is assumed that the material is homogeneous, so that 
material properties are isotropic. 

Figure 1 defines the resultant forces and moments per unit length on an element of a spherical shell under an 
axisymmetric load. Because the shell thermal loading is assumed to be axisymmetric, the change in temperature is 
only a function of (p and z. The area of the element in figure 1 is, therefore, R 2 d(p-d0. The radius of curvature, R, is 
to the middle surface. The resultant forces and moments per unit length are defined below as 


h/2 


h/2 


h/2 


N (p= I ° w - dz N e= l°00- dz M (p=- I z -° w - dz M 0 = 

-h/2 -h/2 -h/2 


h/2 

- jz-o 00 -dz 
-h/2 


h/2 

Q = - f o • dz 

9 -i/2 9Z 


The equilibrium equations for this spherical shell element are given by 


(1) 


d_ 

dip 

d_ 

Stp 


N cp-sin(q>) — N 0 -cos((p)-Q(p -sin((p) = 0 


Qm-sinfa) +N<p -sin((p)+N 0 •sin(cp) = 0 


(2) 

( 3 ) 
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d_ 

ftp 


M e •cos((p)-Q, p -R - sin(tp) = 0 . 


(4) 



•sin((p) 


A more complete derivation of these equilibrium equations can be found in references 1 1 and 12. 

As shown in figure 1, the displacement v is in the tangential direction of increasing (p and the displacement w is 
normal to the surface (i.e., it coincides with the z axis, or radial axis, in a spherical coordinate system). The 
deformations s w and 899 are given by 


(5) 


where 


£ W T, 


dv | 0 1 r l 

— + w £ O0 = — • |y • cot(tp) + wj 
dip ) R 


K m = — 

9 R 


dip 


Kfl = 


X • cot(tp) 
R 


( 6 ) 


and 



dw ' 

v . 

I d( pJ 


(7) 


Here, and Egg are the deformations of the spherical shell’s middle surface. 

The aforementioned deformations and stresses are related by Hooke’s law. The deformations as a function of 
stresses and are given by 


- — • (a tptp - v • CTge)+ a • T Egg — — • (a 00 - v • a (p(p )+ a • T . 


( 8 ) 


Here, T is defined as the change in temperature and is a function of tp and z. To simplify this analysis, it will be 
assumed the change in temperature is linear with thickness so that 


T(<p, z) = i 0 (<p) + z • i(<p) . 


(9) 


Throughout the remainder of this paper, the change in temperature will be referred to as the temperature. 

The deformation equations of eq. (8) can be solved for the stresses. These stresses can then be used to solve for 
the resultant forces and moments per unit length in eq. (1). Equation (7) can be recast with eqs. (5) and (6) to yield x 
as a function of the middle surface deformations. By manipulating these equations further and introducing a new 
variable defined below 


q = 4-Q (p • — , 

v h 2 


(10) 


the following two second order differential equations can be derived for a spherical shell of constant thickness under 
an axisymmetric thermal load 


+ cot(cp) ■ - cot 2 (<p) • X ~ v • x + 

dtp" dtp 


h 2 -R 


4-N 


d^q , dq ? . , 

— - + cot(tp) • cot" (tp) • q + v • q 

dtp" dtp 


4 • E • R 
h 


= (l + v)-a-R- — 
dtp 

(11) 

4 • E • R dr 

• X = a . 

h dtp 

(12) 
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Here, N is defined as 


N = 


E-h 3 

12 • (1 - v 2 ) 


(13) 


The general solutions for x and q are the sum of their non-homogeneous and homogeneous solutions to the second 
order differential equations of eqs. (11) and (12) 

x = x + x q = q + q- (14) 


The non-homogeneous solutions will be solved first. To simplify the analysis, it will be assumed that the change 
in temperature through the thickness as a function of (p is constant. That is, although the temperature can vary 
through the thickness, the temperature variation as a function of thickness is the same everywhere on the spherical 
dome. As a result, di/dcp = 0, so that the right hand side of eq. (11) is zero. To solve these differential equations, the 
first derivative of the function x 0 is expanded in a series with respect to the Legendre, or spherical, functions of the 

first kind, P^fcosfq))] (ref. 13), so that 


dt 


- = S a " -P^osM]. 


dtp , 
x n=l 

The non-homogeneous solutions for eqs. (11) and (12) are, therefore 


X 


= 2 • a T“' A n •Pn 1) [cOS(tp)] 


n=l 


where 


q = 2- a ~ ^B n • PHcos(cp)] 


n=l 


A„ =-- 


1 + V 


{Ml-4f »„} 


B n = ~T~ { [i _ v - 11 • ( n + 1)] • a n } 

A 


A = [l - n • (n + l)] 2 - v 2 + 12 • (l - v 2 )• 


_ 2 ■ n + 1 (n- 1 ) f dT o 


( n { • f -p- ' sin(tp) • Pp [cos(tp)] • dtp . 
(n + 1)! J dtp 


(15) 

(16) 

(17) 

(18) 

(19) 

(20) 
(21) 


The homogeneous solutions of eqs. (11) and (12) are slowly convergent hypergeometric series (refs. 9, 11, 
and 12). For very thin shells where tp e is not small, however, boundary loadings only affect stresses in the vicinity of 
the edge of the spherical dome (refs. 9 and 11). These stresses decrease rapidly toward the center of the shell. 
Therefore, it is possible to neglect y and dx/dtp when compared to d 2 >pdtp 2 in eq. (11), and q and dq/dtp when 
compared to d 2 q/dtp 2 in eq. (12) (refs. 9, 11, and 12). With these assumptions, the homogeneous solutions to 
eqs. (11) and (12) for a closed spherical shell become 
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X = C • e k '^ e ^ • cos[k • (tp e - tp) + ip] 


(22) 


q 


4 • N d 2 x 
h 2 • R dtp 2 


N -, C k ■ e k '^ e ^ -sin[k-(tp e -tp)+ip] 


(23) 


where 


k = 


3 • (1 - v 2 ) • 



(24) 


The constants C and \p are determined by applying an appropriate set of boundary conditions, which will be 
discussed later. 

With the general solutions for q and x above, the stresses per unit length can be solved. The solution for Q^, is 
easily obtained from eq. (10). From the equilibrium eqs. (2) and (3), N,, and N fl can be solved with Q,p 


N (p = “Q<p • cot ( ( p) = 


h 2 


4 • R 


• q • cot((p) 


(25) 


N d Q<P _ h 2 dq 
® dtp 4 • R dtp 


Note that dq/dtp is merely 


= ^q + ^q 

dtp dtp dtp 


(26) 


(27) 


where 


and 


dq , R 

— = 2-a- — 
dtp h 


• X B n • l n ' ( n + J )' P n [cOS(tp)]-COt(tp)‘ P,$ 1) [cOS(tp)]} 

n=l 


(28) 


~~ ~ ~ 8 N c k -e k '^ e ^ -{sin[k-(tp e -tp) + \)/]-cos[k-(tp e -tp)+i|/]}. 
dtp Rh 2 

The moments are given by 

M tp = -N • [k^ + v • kq - (l + v) - a • t] = - ^ • 

M 0 = — N - [klq + v - ic^p - (l + v)- a • T ]- • 

K. 

Further note that d^/dcp is merely 


co 


t(tp)-x + v ~-R-(i + v)- 
dcp 


a • x 


— + v • cot(cp) • x ~ R • (l + v) • a • x 
dcp 


(29) 


(30) 


(31) 
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where 


( 32 ) 


d X _&X t dx 
dip dip dip 


• {n • fn + l)- P n [cos(cp)]— cot((p)* [cos(cp)]} (33) 

and 

-^ L = C-k-e -k '( iPe-ip ) -{cos[k-((p e — cp)+ \|/] + sin[k - (tp e -(p)+i|/]}. (34) 

dcp 

The average stresses o w , Oee, and o< pz given in eq. (1) are merely their respective forces per unit length divided by 
the thickness, h. 

III. Analytical Model of Displacements of a Spherical Shell with an Expanding Fixed Edge 

This section presents how the constants C and \|/ for the above analysis are determined by applying an 
appropriate set of boundary conditions. Although Nowacki determined analytical solutions for the above stresses, he 
did not solve for the tangential and normal displacements (i.e., v and w, respectively). He did, however, provide a 
process for solving the displacements (ref. 9). The following presents this process and the resulting analytical 
solutions. Afterwards, a simplified analytical solution for the displacement at the center of the dome, w 0 , is 
provided. Finally, the process for determining stresses and displacements using the analytical solutions of this study 
are presented. 

A. Boundary Condition and Determination of Constants 

The homogeneous equations for y_ and q can now be solved for a specific set of boundary conditions. While 
these equations can be solved for any number of boundary conditions, a fixed (or clamped) edge support was used 
for this analysis. With a fixed edge support, spherical dome edge rotation and displacements in the direction of the 
chord radius at the edge of the spherical dome are prevented. To make the analysis results more applicable, it was 
further assumed that the edge support at r 0 could be thermally expanded by Ar 0 as shown in figure 2. During this 
expansion, the radius of curvature at the edge support also expands accordingly. Such an assumption is appropriate 
for spherically domed grids that have a fixed dome support that can expand from thermal heating. While such an 
assumption appropriately accounts for the thermal expansion of the grid mounting structure, this edge support is 
assumed to be infinitely rigid, which is not typically the case with a true ion optics mounting system. Unfortunately, 
mounting system designs vary, and the specific mounting system design will affect the boundary conditions that 
should be used in this analysis. The results of this analysis will, therefore, be specific to the selected boundary 
conditions above. However, the following analysis can easily be repeated for a new set of boundary conditions. 

Regarding the temperature, a further assumption made here is that the temperature is constant throughout the 
thickness of the grid. Given that most grid thicknesses are small relative to the radius of curvature and that the 
material is perforated, the assumption is reasonable. So, the temperature distribution is given by 

T(<p, z ) * T(q>) = t q (<p) . (35) 

It will be shown that the boundary conditions for thermally expanding fixed edge support are %(<p e ) = 0 and 
See(<Pe) = Ar 0 /r 0 . With the above boundary condition, the constants C and \|/ in the homogeneous solutions of eqs. (22) 
and (23) can be determined. At <p = cp e , it can be shown that 


v(tp e ) = Ar 0 • cos(cp e ) 

(36) 

w((Pe) = Ar o ’ sin((p e ) 

(37) 


— = 2-a- — -i^Aj, 
dip h ^ 


n=l 
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(38) 


dw(y e ) 

dcp 


= Ar 0 -cos(cp e ) 


Using the results of eqs. (36) and (38) in eq. (7), it is found that x(<Pe) = 0. When the non-homogeneous and 
homogeneous solutions of eqs. (16) and (22), respectively, are used to solve % at cp = <p e , the constant C is found to be 


jCOU 

cos(i|/) 


(39) 


Because the edge support can expand by Ar 0 , the strain at the edge support is, by definition, £g e ((p e ) = Ar 0 /r 0 . 
Equations (1), (25), and (26) can be used in eq. (8) of £99 to yield 


— = £ee(<P e ) = ^ 

i' E 


"'k.v. V 

v h h J 


+ a •!„=-- 


E • 4 • R 


^ - q • v • cot(cp) | + a • t 0 . 

dip 


(40) 


The non-homogeneous and homogeneous solutions of eqs. (17) and (23), respectively, are used to solve for q and 
this, along with eq. (36) for v((p e ), can be used in eq. (40) above. The constant \| / is, thus, found to be 


\|/ = arctan 


R J 


f dq(q> e ) 


4 • R ^ dtp 


— q(tP e ) ' V ' COt(tp e ) +h-E 


Ar„ 


a-x 0 ((p e ) ° 


+ 2-N-k 3 -x(9 e ) 


2 ' N • k 2 ■ [k — v • cot(tp e )] • x(*Pe ) 


(41) 


Note that the constants C and \\i above are for a fixed edge support that can expand by Ar 0 . To obtain a solution for a 
fixed edge with no expansion, the edge expansion term Ar 0 is merely set to zero. 

B. Derivation of Tangential and Normal Displacements 

With the aforementioned equations for forces per unit length solved, the displacement of the spherical dome in 

the v and w directions, shown in figure 1, can now be solved. The middle surface strain equations for e® and 
See of eq. (6) can be combined to yield 


— -v-cot(cp) = R-(£® (p -£e 0 ) . ( 42 ) 

The strain equations e^, and £99 of eq. (8) can be recast to give the stresses and 099 as a function of the strains £ 9(p 
and £99. These stresses can be integrated in the first two integrals of eq. (1) to yield N 9 and N9 as a function of the 

middle surface strains e® and £00 (ref. 9) 

N (P =-^ L -T-[ j 4<p + v ' £ ee -(l + v)-a-T 0 ] (43) 

Y 1 - v 

N 0 =^y'[see + v - £ w -( 1 + v )' a ' T o] • (44) 

Equations (43) and (44) above can be recast to give e^ and £90 as functions of and N e . These can then be 
used in eq. (42), along with eqs. (25) and (26) for and N e , respectively, to yield 
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dv ./ , R • (1 + v) (1 + v) • h 


dq , , 

— - q • cot(tp) 
dtp 


The differential equation above can be solved by assuming a general solution 

(1 + v) • h 


v(<p) = - 


4 • E 


■ q((p)+ f(<p) 


(46) 


where f(tp) represents the homogeneous solution portion of the general solution. The general solution above can be 
used in eq. (42) to yield 


f • cot(cp) = 0 

dcp 


whose solution is 


f(cp) = C v • sin(cp) . 


(47) 


(48) 


The general solution for v is, therefore 


(1 + v) • h 
4 • E 


•q + C v • sin(cp) . 


(49) 


The solution for the displacement w can be found substituting eq. (49) for v back into the middle surface strain 
equation for Se 0 of eq. (6). Using the equation for Egg as a function of N 9 and N 0 , along with eqs. (25) and (26) for 
N 9 and N e , respectively, it can be shown that 


w = Rai„ -- 


4 • E 


dq . , 

— + q-cot(tp) 
dtp 


-C v -cos(tp) . 


(50) 


The constant C v can be determined by noting that at cp = tp e , v is given by equation (36). It can be shown that 


C v = 


Ar 0 • cos(tp e ) - (l + ^l h • q(9e ) 
4E 

sin(tPe) 


(51) 


C. Simplified Analytical Expression for the Displacement w 0 at the Spherical Dome Center 

The displacement w at the center of the spherical dome, or tp = 0, is of interest because this location experiences 
the largest displacement for a temperature distribution that monotonically decreases with increasing <p. The 
displacement can be simplified by noting that as cp approaches zero, q in equation (50) for w goes to zero. This is 

because the non-homogeneous solution of q in eq. (17) is zero since P d) (1) =0 and the homogeneous solution of q in 
eq. (23) is negligibly small because it is multiplied by cxpf— k tp e ), which is a small value for large ratios of radius of 
curvature to thickness. Furthermore, it can be shown that dq/d(p also becomes negligibly small as (p approaches zero. 
This is because the homogeneous solution of dq/dcp in eq. (29) is also negligibly small since it is also multiplied by 
exp(-k(p e ) and the non-homogeneous solution of dq/d(p in eq. (28) typically produces a negligibly small value near 
the center. 

These findings are not too surprising. The terms q and dq/dcp are directly proportional to o w and Gee, 
respectively. One benefit of spherically-domed shells under thermal loads is that these stresses tend to decrease 
significantly near the center of the dome, as opposed to flat plates, whose stresses tend to increase under similar 
thermal loads. This is because the spherical domes have a preferential direction in which to thermally expand. 

So, the resulting displacement w at (p = 0 reduces to 
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w 0 « R • a • T(o)- C v . (52) 

The first term on the right hand side of the above equation represents the thermal expansion of a spherical shell of 
radius R, thermal expansion coefficient a, and temperature T at <p = 0. The second term (i.e., C v ) represents the 
effect of the boundary condition, given in eq. (51). 

The expression for C v is complicated because it includes q(<p e )- It is, however, possible to simplify this 
expression. The homogeneous solution to q is given by combining eqs. (23) and (39) at q> = q> e to give 


q(<Pe) = 


8-N-k^ 
h 2 R 


■ tan(\[/) • %(<Pe ) 


Using this, along with eq. (41) for \| i, the general solution to q(cp e ) becomes 


q(tPe) = q( ( Pe) + - 


dq((Pe) - 

dip 


4RE 


Ar n 8-N-k J 


Q ) , , , i XV i—i i-u 0 O i >1X 

q((p e )-v-cot((p e ) + — a-T 0 ((p e ) + x(*Pe ) 

n r T) U ^ 


k-v-cot(cp e ) 


Using eqs. (16), (17), and (28) in the above equation, it can be shown that 


q(tp e ) = - 


4-R-E f Ar 0 

h l 

k-v-cot((p e ) 


q oo loci 

2 ' a ' TT ' Z! ] Bn ' F + ))• p n !) [cos((p e )]- n • (n + 1)- P n [cos(cp e )]J+ " • A n • [cos((p e )] 


•11- V 2 , 


k-v-cot((p e ) 


For thin shells, R/h is large when compared to n 2 in the series summation, so eq. (20) for A reduces to 


(i-vd 


A « 12-11- v • — 


so that eqs. (18) and (19) for A n and B n , respectively, reduce to 


A « a 

2 • R 


[l - v - n • (n + 1)] • 

4-4fAf 


These equations can now be used in eq. (55). It can be shown that the first term in the series summation of eq. (55) is 
always negligibly small compared to the second term. Equation (55), therefore, reduces to 


4 • R • E [ , , Ar 0 2-a-E-k ^ n(1)r , q 

a-T 0 (9e) , / f / .an'Pn F S (<Pe) 

h l O h- i-v 2 t, 


k - v • cot((p e ) 
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But the series summation of the above equation is merely the temperature gradient of eq. (15). Using this with the 
temperature assumption of eq. (35) in the above equation yields 


4 • R • E 
h~ 


a-T(cp e )- — 


q(9e) ! 


2 • a • E • k dT(cp e ) 

M-v 2 )' d(p 


k-v-cot(cp e ) 

With this simplified expression for q at (p = <p e , eq. (60) can be used in eq. (51) for C v to yield 


(1 + v)- R 


a-T(cpe)- — 


q-k-h dT(tp e ) | 


C v ~ Ar 0 -cot((p e )-- 


■yJl2-( l-v 2 ) d( P 


k • sin(tp e ) - v • cos((p e ) 


(60) 


(61) 


Equations (52) and (61) define a simplified expression for the displacement of the center of the spherical dome 
requiring knowledge of material properties, geometry of the dome, temperatures at the center and edge, and the 
derivative of the temperature at the edge. As expected, the deflection of the dome center is a function of the 
material’s Poisson’s ratio but is independent the material elastic modulus. 

D. Solution Process for Thermally-Induced Moments, Stresses, and Displacements 

The process for solving for the thermally-induced displacements of the spherically domed shell, along with the 
stresses and moments, is presented below. The required inputs include shell material properties, geometry, 
and temperature distribution. Shell material properties include a, v, and E. Shell geometry inputs include R, h, and 
(p e . Often, however, only the dome height and chord radius are known. The radius of curvature can then be 
determined by 


R = 


r 2 +h 2 
i o ' 11 dome 

2 ’ hdome 


while (pe is determined by 


/ 

(p e = arcsin • 

v 


R 


A 


J 


(62) 


(63) 


The temperature distribution must be some continuous mathematical function of (p. Finally, a value for Ar u 
is required. 

The derivative of the temperature distribution, dT/d(p = di./dcp. is first determined. This result is used in eq. (21) 
to determine the values for a n . Solutions to the spherical functions of the first kind, along with the Legendre 
functions, can be defined using known Legendre explicit expressions and recurrence relations for higher orders and 
degrees (ref. 13) The results for a n , along with A from eq. (20), are used in eqs. (18) and (19) to determine values for 
A n and B n , respectively. These results are then used in eqs. (16), (17), (33), and (28) to determine the non- 
homogeneous solutions of %, q, dx/d<p, and dq/d<p, respectively. The results of these four equations are used with the 
temperature at (p = cp e to determine the constants \p and C of eqs. (41) and (39), respectively. With these constants 
and eqs. (24) and (13) for k and N, respectively, the homogenous solutions to %, q, dx/d<p, and dq/d(p, in eqs. (22), 
(23), (34), and (29), respectively, can be determined. The general solutions of x, q, dx/dtp, and dq/d<p are, therefore, 
the sum of their respective non-homogeneous and homogeneous solutions. With these solutions, Q, p . N„ and N, :) are 
determined with eqs. (10), (25), and (26), respectively. The moments M,, and M 0 are determined with eqs. (30) and 
(31). The average stresses o w , o ee , and o 9Z are merely the forces per unit length divided by the thickness, h. 

Regarding displacements, the constant C v in eq. (51) is determined with the general solution of q. With this 
result and the general solution for q, the tangential displacement, v, is solved with eq. (49). The general solutions of 
q and dq/d<p of are used to determine the normal displacement, w, in eq. (50). 
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There are limitations to the solutions to N e , N ip , M<p, M 0 , and w. Specifically, the solutions to these equations 
produce a singularity at 9 = 0, which is the geometric center of the spherical dome. This is due to the cot ( 9 ) term in 
each equation. That the cot(ip) term was not properly accounted for in the solutions is an artifact of the assumption 
made for the homogeneous solutions of % and q. However, solutions can be obtained in the immediate vicinity of 
(p = 0 and the results extrapolated to the dome center. 

The non-homogeneous solutions of x and q in eqs. (16) and (17), respectively, and their derivatives include a 
series summation to infinity. These summations converge after some number of summations, so that some limiting 
value other than infinity can be used. To determine this value, it is important to understand that the non- 
homogeneous solutions to the differential equations of eqs. ( 11 ) and ( 12 ) were found by assuming that the derivative 
of the temperature distribution across the spherical dome could be expressed by the series summation in eq. (15). 
The required number of summations is, therefore, the number of summations necessary for eq. (15) to accurately 
express the actual derivative of the temperature as a function of 9 . For larger values of n, a n becomes negligibly 
small to ensure a convergent solution. So the required number of summations can be determined by comparing the 
actual derivative of the temperature distribution with that of eq. (15). 

IV. Effects of Shell Perforation 

Because ion optics are perforated, the effect of this perforation on the displacements must be included in the 
analytical model. Fortunately, determining the stresses and strains of perforated plates has been dealt with 
considerably in literature, examples of which are included in references 6 , 14 to 17. The effects of plate perforation 
on displacements and stresses are determined by modeling the perforated plate as an equivalent solid plate with 
modified, or effective, material properties. The effective material properties include the elastic modulus and the 
Poisson’s ratio (ref. 15). The resulting effective material properties are functions of the perforation pattern layout, 
aperture center-to-center spacing, and aperture diameter (ref. 14). With these effective material properties, 
displacements and nominal stresses are determined. The displacements in the actual perforated plate are the same as 
those of the equivalent solid plate utilizing effective material properties. Actual ligament stresses are determined by 
multiplying the nominal stresses with predetermined stress multipliers (ref. 14). 

While this process offers significant simplicity, a short-coming is that effective material properties can be a 
function of the type of load applied to the perforated plate. This is especially true for thinner plates with a high 
perforated open area fraction and where the ratio of the material thickness divided by the center-to-center aperture 
spacing is small (ref. 15). The resulting equivalent solid plate possesses effective material properties that are not 
isotropic, which violates a basic assumption used in developing theoiy of thermal displacements in thin-walled 
shells. A thermally-loaded spherically domed grid will have both in plane stresses (i.e., thermally-induced 
compressive and tensile loads in the plane of the perforation), and bending loads. Therefore, if effective elastic 
moduli and Poisson’s ratios vary with the type of load, this will introduce errors in determining displacements. 

Effective material properties that account for perforation can be determined theoretically (refs. 14 to 17) or via 
finite element modeling (ref. 6 ). The latter technique is labor-intensive and can be shown to produce results that are 
similar to those of the theoretical study of reference 14. The theoretical studies of references 14 to 17 demonstrated 
that the ratios of effective to actual E and v are functions of the minimum webbing width divided by the center-to- 
center aperture spacing for a given aperture layout. These studies also provided plots so that effective material 
properties could be readily determined for a given ratio of minimum webbing width to center-to-center 
aperture spacing. 

V. Comparison of Analytical Model with Finite Element Model Results 

The analytical model above was compared to the results from the finite element model in reference 6 by Shunk. 
Comparisons were made to both solid and perforated spherical shells with varying geometric and material 
properties. This was done so that the efficacy of the analytical solution and the process for accounting for 
perforation could be investigated separately over a relatively broad range of material and geometric properties. 
Furthermore, comparisons were only made for displacements normal to the surface because only these data were 
available for comparisons. Regardless, some tangential displacements from the analytical model will be presented. 

Determining the efficacy of this analytical model by comparing it to a finite element model was done for three 
reasons. First, comparisons were direct because inputs such as boundary conditions and temperature distributions 
could be exactly matched (i.e., both used fixed edge supports). Second, comparisons could be made throughout the 
spherical dome because these data could be easily determined from the finite element model. Finally, measured data 
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with similar boundary conditions and free of other factors that could obfuscate comparisons could not be found by 
this author. 

A. Solid Shell Normal Displacements 

The finite element model in reference 6 was used to calculate the thermally-induced displacements normal 
to the spherical surface, w, as a function of chord radius for three different materials (i.e., 304 stainless steel, 
titanium, and molybdenum) and for two different chord radii (i.e., 50 and 75 cm). Although the aforementioned 
analytical model yields results as a function of <p, this variable can be converted to a chord radial dimension with the 
following relation 


r = R-sin(tp). (64) 

For this comparison, the spherical domes were assumed to be solid, and not perforated. Material properties are listed 
in table 1. The screen grid thicknesses of table 2 were used. Geometric properties of the screen grid spherical domes 
listed in table 3 were used. The temperature distribution used here was that developed in reference 6. The 
temperature as a function of chord radius was a second order quadratic expression 

T ((p) = Ajq + Ajj R- sin(cp) + A j 2 * R ~ * sin ~ ((p) . (65) 

The constants used in the above temperature equation are those for the screen grids listed in table 4. Finally, the 
boundary condition in reference 6 was a clamped edge with no chord expansion, so that Ar 0 was equal to zero in the 
analytical model. 

The normal displacements from the analytical and finite element models are shown in figures 3 and 4 for the 50 
and 75 cm chord radii, respectively. As the figures show, there is excellent agreement between the finite element 
model of reference 6 and the analytical model of this study. Table 5 compares the simplified analytical solution of 
eqs. (52) and (61) for w G to the finite element model results from reference 6 at the center of the solid spherical 
dome. The comparisons were excellent, with the worst-case difference being < 4 percent. 

B. Perforated Shell Normal and Tangential Displacements 

For the perforated material, only the 50 cm chord diameter was examined. Two aperture diameters were 
examined, and these were labeled screen and accelerator grid. The screen and accelerator grid aperture diameters 
were 0.899 and 0.399 cm, respectively, with an aperture center-to-center spacing of 0.953 cm for both. The same 
three materials as the prior section were examined. Effective material properties due to shell perforation under a 
plane stress load were determined from reference 14 and are listed in table 6. The grid thicknesses of table 2 were 
used. Geometric properties of the spherical domes are listed in table 3 for the 50 cm chord diameter. The constants 
used in the temperature equation are listed in table 4 for the 50 cm screen and accelerator grids. The equivalent finite 
element model results for the perforated grids were taken from Shunk, who modeled a spherical shell with apertures 
(ref. 6). Because of symmetry, Shunk modeled a perforated quarter section of the domed spherical grid using 
appropriate boundary conditions. 

The normal displacements from the analytical and finite element models are shown in figures 5 and 6 for the 
screen and accelerator grids, respectively. Table 7 compares the simplified analytical solution of eqs. (52) and (61) 
for w 0 to the finite element model results from reference 6 at the center of the perforated dome. As the figures and 
table show, there was good agreement between the analytical and finite element models for the titanium and 
molybdenum materials. The stainless steel material produced poorer comparisons, however. 

Tangential displacements using the analytical model were determined for the perforated molybdenum and 
titanium grids. These tangential displacements are shown in figure 7 for both grids. Tangential displacements at the 
center and edge were both zero due to symmetry and setting Ar 0 = 0, respectively. 

C. Discussion of Results 

For the solid shell results, normal displacement comparisons showed that the analytical model produced results 
that closely matched the finite element model of reference 6. This is not surprising because the assumptions used 
developing the above equations for thermally-induced shell displacements were the same as those commonly used 
for thin-walled shells of revolution (ref. 11). The simplified equation for the normal displacement of the spherical 
dome center was also found to accurately predict the displacement in this location, demonstrating that the simplified 
equation for w Q can be an accurate predictive tool. 
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Perforated shell normal displacement comparisons demonstrated that the use of effective material properties can 
produce inaccurate results in some cases, under-predicting the displacements at and near the grid center. It is 
speculated the discrepancies were likely the result of the steel’s material properties and the use of effective material 
properties that can vary with load type. Reference 15 showed that a perforated plate with ratios of the material 
thickness divided by the center-to-center aperture spacing similar to all of those of this study should produce 
effective material properties that are a function of load type. It can be shown that the moments of eqs. (30) and (31) 
are directly proportional to 


.. a-E-h 

oc ocMa. (66) 

9 R 9 

Because these moments produce bending loads, the steel screen grid would experience the highest bending loads 
because it had the highest thermal expansion coefficient, the second highest elastic modulus, and the highest 
thicknesses. As a result, the stainless steel produced poor comparisons to the finite element model. 

The analytical model results for normal displacements still produced accurate results for the titanium and 
molybdenum materials, likely due to their lower thermal expansion coefficients. Table 7 shows that the simplified 
equation for w G yielded results that were within 3 percent of the results from the finite element model. Because most 
grid materials are either molybdenum or some other material with a low thermal expansion, this analytical approach 
should produce accurate results for these materials. 

The analytical model results for tangential displacements showed that tangential displacements were 
significantly smaller than normal displacements. While it is likely the case that normal displacements will always be 
greater that tangential movement, the significantly small values in figure 7 were due, in part, to setting Ar 0 = 0. A 
larger value for Ar 0 would have caused larger tangential displacements. 

From the aforementioned perforated plate results, changes in grid gap and grid-to-grid aperture alignment can 
be easily determined. At a given location r, grid gap and grid-to-grid aperture alignment changes are merely the 
difference of the normal and tangential displacements, respectively, of each grid. 

VI. Conclusions 

An analytical method for calculating the thermally-induced normal and tangential displacements of spherically 
domed grids was described. The theory for thin-walled shells of revolution under thermal loads varying both through 
the thickness and along surface was reviewed. Only spherical shells of uniform thickness and under an axisymmetric 
thermal loading were considered. A fixed edge support was used for this analysis. To make the analysis results more 
applicable, it was further assumed that the edge support could be thermally expanded. 

Displacements both normal and tangential to the surface of the spherical shell were derived. A simplified 
equation for the displacement at the center of the spherical dome was also derived. This simplified expression for the 
displacement only requires knowledge of material properties, geometry of the spherical dome, temperatures at the 
center and edge, and the derivative of the temperature at the edge of the dome. 

The effects of plate perforation on displacements and stresses were determined by modeling the perforated plate 
as an equivalent solid plate with modified, or effective, material properties. The effective material properties 
included the elastic modulus and Poisson’s ratio. The resulting effective material properties are functions of the 
perforation pattern layout, aperture center-to-center spacing, and aperture diameter. The displacements in the actual 
perforated plate are the same as those of the equivalent solid plate utilizing effective material properties. 

The analytical model was compared to the results from a finite element model for displacements normal to the 
surface. For the solid shell results, comparisons showed that the analytical model produced results that closely 
matched the finite element model results. The simplified equation for the normal displacement of the spherical dome 
center was also found to accurately predict the displacement in this location. For perforated shells, the analytical 
solution and simplified equation produced accurate results for the titanium and molybdenum materials, likely due to 
their lower thermal expansion coefficients. Because most grid materials are either molybdenum or some other 
material with a low thermal expansion, this analytical approach should produce accurate results for these materials. 

The results from this study demonstrated that this analytical model can accurately predict grid displacements 
and can, therefore, be used as a predictive tool. The analytical equations presented in this study are further useful 
because they provide the relevant geometric and material parameters and their impact on grid displacements. 
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TABLE 1.— SPHERICAL DOME MATERIAL 


PROPERTIES AND THICKNESSES FROM REF. 6 


Material 

E, 

V 

a, 

GPa 


m/m-°C 

304 stainless steel 

193 

0.3 

17.8xl0“ 6 

Titanium 

100 

0.33 

9.5xl0“ 6 

Molybdenum 

300 

0.3 

5.43xl0 -6 


TABLE 2.— SPHERICAL DOME THICKNESSES FROM REF. 6 


Material 

Screen h, 
cm 

Accelerator h, 
cm 

304 stainless steel 

0.300 

0.396 

Titanium 

0.152 

0.203 

Molybdenum 

0.152 

0.396 


TABLE 3.— SPHERICAL DOME GEOMETRIES FROM REF. 6 


Chord diameter 

R, 

cm 

<Pe 

50 cm 

85.75 

16.95° 

75 cm 

128.62 

16.95° 


TABLE 4.— TEMPERATURE DISTRIBUTION CONSTANTS FROM REF. 6 


Chord diameter and grid 

Ato 

An 

A T 2 

50 cm screen 

-0.095 

-0.636 

283.3 

50 cm accelerator 

-0.096 

-0.356 

237.3 

75 cm screen 

-0.042 

-0.423 

283.4 


TABLE 5.— COMPARISON OF SIMPLIFIED SOLUTION w 0 TO SOLID 


SHELL FINITE ELEMENT MODEL RESULTS FROM REF. 6 


Material 

50 cm chord radius 

75 cm chord radius j 

FEM, 

cm 

Analytical, 

cm 

Error, 

% 

FEM, 

cm 

Analytical, 

cm 

Error, 

% 

304 stainless steel 

0.504 

0.504 

0.2 

0.723 

0.735 

1.7 

Titanium 

0.259 

0.258 

-0.3 

0.376 

0.379 

0.9 

Molybdenum 

0.145 

0.147 

1.6 

0.210 

0.216 

3.2 


TABLE 6.— PERFORATED SPHERICAL DOME EFFECTIVE MATERIAL PROPERTIES 



Screen grid 

Accelerator grid | 

Material 

E, 

V 

E, 

V 


GPa 


GPa 


304 stainless steel 

5.02 

0.795 

124 

0.315 

Titanium 

2.60 

0.795 

64.1 

0.334 

Molybdenum 

7.80 

0.795 

192 

0.315 


TABLE 7.— COMPARISON OF SIMPLIFIED SOLUTION w 0 TO PERFORATED 


SHELL FINITE ELEMENT MODEL RESULTS FROM REF. 6 


Material 

Screen grid 

Accelerator grid ! 

FEM, 

cm 

Analytical, 

cm 

Error, 

% 

FEM, 

cm 

Analytical, 

cm 

Error, 

% 

304 stainless steel 

0.655 

0.573 

-12.5 

0.447 

0.432 

-3.47 

Titanium 

0.289 

0.281 

-2.77 

0.221 

0.219 

-0.90 

Molybdenum 

0.165 

0.161 

-2.42 

0.136 

0.132 

-2.94 
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Figure 1. — Element of a spherical shell with the resultant forces and moments per unit length 
from an axisymmetric thermal load. Radius of curvature and resultant forces 
and moments are to the element middle surface. 



Figure 2. — Illustration of the boundary conditions used in this study for a spherically domed grid. 
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Figure 3. — Solid spherical dome displacements normal to the surface as a function of chord 
radius for the 50 cm chord diameter dome. Finite element model (FEM) 
results are from reference 6. 
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Figure 4. — Solid spherical dome displacements normal to the surface as a function of chord 
radius for the 75 cm chord diameter dome. Finite element model results are from reference 6. 
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Figure 6. — Perforated spherical dome displacements normal to the surface as a function 
of chord radius for the 50 cm chord diameter accelerator grid. Finite element 
model results are from reference 6. 
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Figure 7. — Perforated spherical dome displacements tangential to the surface as a 
function of chord radius for the 50 cm chord diameter screen and accelerator grids. 
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results from a finite element model. For the solid shell, comparisons showed that the analytical model produces results 
that closely match the finite element model results. The simplified equation for the normal displacement of the spherical 
dome center is also found to accurately predict this displacement. For the perforated shells, the analytical solution and 
simplified equation produce accurate results for materials with low thermal expansion coefficients. 
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